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Abstract 



We describe a measurement of the charge asymmetry of leptons from W 
boson decays in the rapidity range < \yi\ < 2.5 using W— > eu, \iv events 
from 110±7 pb" 1 of data collected by the CDF detector during 1992-95. The 
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asymmetry data constrain the ratio of d and u quark momentum distributions 
in the proton over the x range of 0.006 to 0.34 at Q 2 ~ M^. The asymmetry 
predictions that use parton distribution functions obtained from previously 
published CDF data in the central rapidity region (0.0 < \yi\ < 1.1) do not 
agree with the new data in the large rapidity region (\yi\ > 1.1). 



PACS numbers: 13.85.Qk, 13.38.Be, 14.70.Fm 
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At Tevatron energies, W + {W~) bosons are produced in pp collisions primarily by the 
annihilation of u (d) quarks in the proton and d (u) quarks from the antiproton. Because u 
quarks carry on average more momentum than d quarks the W + bosons tend to follow 

the direction of the incoming proton and the W~ bosons that of the antiproton. The charge 
asymmetry in the production of W bosons is related to the u and d quark distributions at 
Q 2 « M 2 ,, and is roughly proportional @|| to the ratio of the difference to the sum of the 
quantities -(xi) and -(^2), where x\ and %i are the fractions of the nucleon's momentum 
carried by the quarks in the p and p, respectively, -(x) is the ratio of the d to u quark 
parton distribution functions for quarks carrying a momentum fraction, x, of the nucleon's 
momentum. 

W asymmetry measurements || from 20 pb _1 of data collected by the Collider Detector 
at Fermilab (CDF) in 1992-93 demonstrated the measurement to be more sensitive than deep 
inelastic scattering experiments to the - ratio within the x range accessible by the W data at 
the Tevatron. These data contribute to global analyses used to extract parton distribution 
functions (PDFs) in the nucleon These functions are used in the calculation of all 

hadronic cross sections. CDF's previous asymmetry data || improved the understanding of 
the u and d quark momentum distributions in the proton in the x range 0.007 < x < 0.24. 
The use of these data resulted in a reduction in the uncertainty from PDFs in the W mass 
measurement |7j from rs 75 to ~ 50 MeV/c 2 . A further reduction of this error to ~ 20 MeV/c 2 
is possible || using the new data presented here. Recently fixed target Drell-Yan data f| 
has also been used to constrain the (u — d) sea quark distribution. 

In this letter, we describe the asymmetry results at y/s = 1.8 TeV using data from the 
entire 1992-95 period, using a total integrated luminosity of 110 ± 7 pb _1 - a 5-fold increase 
relative to the 1992-93 data. The new measurements agree with and supersede the previous 
results. The asymmetry measurement is extended to larger rapidity |TUj (\yi\ < 2.4) by the 



introduction of a new charge determination technique |TT| for electrons, and also by using 



data from the forward muon detector (1.9 < \yi\ < 2.5). These electron and muon data at 
large rapidity provide information about PDFs for a larger x range (0.006 < x < 0.34) than 
previously available. The mean x values of the d and u parton distributions probed by the 
CDF asymmetry data are shown in Table |. 

Since the W rapidity is experimentally undetermined because of the unknown longitudi- 
nal momentum of the neutrino from the W decay, we measure the lepton charge asymmetry 
which is a convolution of the W production charge asymmetry and the V — A asymmetry 
from the W decay. The two asymmetries tend to cancel at large values (\yi\ > 2) of rapidity. 
However, if one assumes W decays proceed via a pure V — A interaction, the lepton asym- 
metry is still sensitive to the parton distribution functions. The lepton charge asymmetry 
is defined as: 

= da + /dyi - da'/dyi ^ 
do + 1 dyi + da~ / dyi 

where da + (da~) is the cross section for W + (W~) decay leptons as a function of lepton 
rapidity, with positive rapidity being defined in the proton beam direction. If the detection 
efficiencies and acceptances for Z + and l~ are equal, then the uncertainties on these quantities 
do not affect A(yi). 
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The tracking detectors at CDF Jl2] are the Vertex Time Projection Chambers (VTX), the 
Central Tracking Chamber (CTC, \rj\ < 1.8), the Silicon Vertex Detector (SVX, \rj\ < 2.3), 
the Central Muon Chambers (CMU and CMX, \r)\ < 1.0), and the Forward Muon Detector 
(FMU, 1.9 < \r)\ < 3.5). The VTX provides r-z tracking information out to a radius of 22 
cm for \r)\ < 3.25 and is used to measure the z location QTUJ of the primary interaction vertex. 
The CTC is an 84 layer drift chamber inside the 1.4 T solenoidal magnet which provides the 
curvature measurement of electrons and muons from W decays. The CTC has a momentum 
resolution of 8(l/pr) = 0.00081 (GeV/c) -1 0. The SVX is a 51 cm long, 4 layer precision 
vertex detector with a track position resolution |TT] of « 10 /im. Because pp interactions at 
CDF are spread along the beamline with an RMS of ~ 30 cm, the geometrical acceptance 
of the SVX is about 60% for W candidate events. In this analysis, the SVX complements 
the CTC for electron tracking at high \r)\ (1.2 < \rj\ < 2.3). The FMU consists of two steel 
magnets (toroids) alternately inserted between 3 drift chambers. 

W decays to leptons are identified by a charged track pointing either at hits in the muon 
chambers or a cluster of energy in the electromagnetic calorimeters (EM) accompanied by 
large missing transverse energy {fir) ||10|| . Electron candidates are required to fall within 
the fiducial regions of either the central (\yi\ < 1.1) or the plug (1.1 < \yi\ < 2.4) EM 
calorimeters and to pass identification cuts based on the EM shower's profile determined 
with test beam electrons. Muon candidates are required to have a track in the muon tracking 
system and a minimum ionizing particle signal in the hadronic and EM calorimeters. The 
transverse energy (Er) of the lepton and Et are required to be greater than 25 GeV. To 
further reduce the background from dijet events where one jet is incorrectly reconstructed as 
a charged lepton and where a large Et results from shower fluctuations or uninstrumented 
regions in the detector, events with a jet [13] with E^ > 30 GeV are rejected. 

The data are divided into five samples: central electrons (\yi\ < 1.1), central muons 
(\yi\ < 1.0), plug electrons within the SVX fiducial region (plug-SVX, 1.1 < \yi\ < 2.4), 
plug electrons outside the SVX geometrical acceptance region but with a CTC track (plug- 
CTC, 1.1 < \yi\ < 1.8), and forward muons (1.9 < \yi\ < 2.5). In total 88047 events are 
selected. A reliable charge determination and a good energy measurement are essential for 
this analysis. Central electrons, central muons, and plug electrons outside the SVX are 
required to have an associated CTC track. A reliable charge determination is ensured by a 
CTC track curvature (C) significance requirement oiC/SC > 2 (where 5C is the error in the 
curvature measurement). The charge of forward muons is measured well in the FMU toroidal 
magnets. The energy of plug electrons is measured by the EM calorimeters, and the rj and 



are measured by the tower segmentation [10|. However, tracking information is needed 



in order to determine the charge. The geometrical acceptance of the plug calorimeters is 
only partially covered by the CTC. Therefore, the efficiency for finding a CTC track, and 
hence determining the charge, drops quickly as a function of rapidity and becomes zero for 
\vi\ > 1-8. 

A new charge determination technique [[ll]] is introduced to identify the charge of plug 
electrons within the SVX geometrical fiducial region. At CDF, positively (negatively) 
charged particle trajectories are bent in increasing (decreasing) <fi inside the solenoidal mag- 
netic field. The SVX track stub measures precisely the initial track direction in 0. The 
plug EM calorimeter (PEM) [|14]] shower centroid measurement combined with the location 
of the vertex position yields another measurement of 0, and the sign of the difference of 
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the two measurements determines the charge. The shower centroid is measured with a 
strip position detector placed at shower maximum in the pseudo-rapidity region between 
1.2 — 1.8. In the pseudo-rapidity region between 1.8 — 2.4, where there is no strip detector 
coverage, the centroid is measured from the calorimeter tower information in r] and 0. We 
define the ratio R = 54> messmed / | ^expected |, where 5(j) mcasuied = PE m - 4>s vx, and 50 exp ected 
is calculated from the calorimeter energy, the radial location of the electron shower and the 
magnetic field strength. R peaks at +1 and -1 for positrons and electrons, respectively. 
The measurement error on R is 0.3, 0.5 and 0.8 for \yi\ bins centered at 1.3, 1.8 and 2.2, 
respectively. This technique effectively doubles the number of plug electrons that can be 
used in the measurement and extends the charge measurement to \y{\ = 2.4. 

A possible charge or Et dependence of each trigger is investigated by using data from 
several independent triggers. No evidence of charge dependence is found in the electron and 
central muon data samples. The forward muon data taken with the toroid polarities set to 
focus n + and respectively, are averaged to cancel out effects of any charge bias in the 
trigger from geometrical acceptance and alignment errors. The efficiency of the plug electron 
triggers does not depend on charge, but is not 100% at Et of 25 GeV. At Et = 25 GeV 
the average trigger efficiency for the plug-SVX sample is 78%. A correction is applied to 
the measured asymmetry in each rapidity bin, using a Monte Carlo calculation and the 
measured trigger efficiency as a function of Et and yi. The correction to A(y{) is found to 
be less than 0.005 compared with a typical statistical error of 0.015 in the PEM region. 

Sources of charge bias in the event selection are investigated by selecting high Et leptons 
from either a sample of Z or W events, which satisfy tight kinematic constraints. No charge- 
dependent bias in acceptance or efficiencies has been found. The charge misidentification 
rate for muons is negligible. The electron charge misidentification rate (e.g. due to multiple 
tracks from conversion of Bremsstrahlung photons) in the central region is determined from 
the rate of same sign central-central Z electrons. In the plug region, the rates are determined 
from W electron data using the number of events outside the two peaks in the R distribution. 
This is checked against the same sign central-plug Z electrons as a function of rapidity. The 
electron charge misidentification rate ranges from 0.2% in the central data sample to 10% in 
the largest rapidity bin {\y{\ = 2.20) of the plug electron data. The charge misidentification 
rate in the FMU sample is < 1%. The effect of the charge misidentification acts to dilute 
the charge asymmetry and is corrected on a bin by bin basis. 

All backgrounds have been investigated and found to be small (see Table |J). The 
average W — > tv background is 2.0 ± 0.2% where the error accounts for the small rapidity 
dependence. In the plug electron sample, the background from misidentified dijet events is 
the largest. This background is charge-symmetric and dilutes the asymmetry. In the central 
electron sample, the background from W — > tv — > evvv is the largest. For the central 
muon sample the largest background is misidentified Z — > where one of the muons is 

outside the acceptance of the CTC. In the latter two cases, the backgrounds are not charge- 
symmetric and the corrections are made using the asymmetry estimated by Monte Carlo 
calculations. The background from misidentified Z decays to electrons is negligible because 
the EM calorimeters have a much larger geometric acceptance than the muon chambers 
or the CTC. In the forward muon sample the background from Z events is 6.5% and the 
charge-symmetric background totals 8% of which half is from dijets and half from lower pt 
muons which are misreconstructed because of random extra hit background in the forward 
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drift chambers. The Z — > t + t~ contamination is negligible. The background from cosmic 
rays in the muon sample is < 0.2%. The A(y{) values are corrected for the backgrounds on 
a bin by bin basis. 

By CP invariance the asymmetry at positive yi is equal in magnitude and opposite in 
sign to that at negative y\ allowing the two measurements to be combined. Figure [I] shows 
the fully corrected asymmetry after taking the weighted mean of the various data sets and 
combining the positive and negative yi bins. The asymmetry and uncertainties are listed in 
Table [iTT[ Figure [3] shows the predictions of Quantum Chromodynamics (QCD) calculated 
to Next-to-Leading-Order using the program DYRAD |T5| with several parameterizations of 
parton distribution functions as input. The MRS-R2 []16| and CTEQ-3M [f2| PDFs have 
been extracted with the inclusion of the 1992-93 asymmetry data in their global fits (the 
1992-93 data constitutes 20% of the sample reported here). As shown in Figure [I], the 
predictions using these PDFs are in good agreement with the present data in the central 
region {\y\\ < 1.1), with x 2 per degree of freedom (d.o.f) values in the range 0.7-1.3. However, 
at high rapidity, these predictions are generally higher than the data, with % 2 /d.o.f values in 
the range 5-10, indicating that the PDF parameterizations should be modified in the range 



0.006 < x < 0.34 (see Table ffl). A recent re-analysis [17] of NMC muon scattering data [18 



on hydrogen and deuterium with improved corrections for nuclear binding effects in the 
deuteron has shown the need for a correction to the - ratio of the form : <5(-) = 0.1(:e + :e 2 ). 
The prediction using the corrected MRS-R2 PDF is shown in Figure p] and shows reasonable 
agreement with the measured asymmetry at high \yi\. Recent global parton distribution 



fits |19| (e.g. MRST) that have incorporated the data presented in this paper also lead to 
agreement at high \y{\. The theoretical uncertainties arising from the effect of the finite 
charm quark mass [^(J are much smaller than the measurement errors. 

The W decay lepton charge asymmetry is mostly sensitive to the - ratio. However, at 
large rapidity it is also affected by the W production spectrum. The DYRAD calculation 
does not reproduce the production pt spectrum of W events at low pt- Therefore, it is 
necessary to check how well the DYRAD calculation compares to one more suitable for the 
low W pt region. The DYRAD prediction is compared to a calculation including soft gluon 



resummation at all orders in perturbation theory implemented in the program REB0S [21 



Both programs yield identical results for the W rapidity distribution, but different results 
for the W pt distribution. At Tevatron energies soft gluon radiation is mainly responsible 
for the p T of W and Z bosons in the range of Pt < 30 GeV/c, and REB0S can reproduce 
the pt spectrum of Z bosons ||22|| . Figure ^ shows the comparison between the asymmetry 



predictions of the DYRAD and REB0S programs using CTEQ-3M |2J PDFs. The two calcu- 
lations agree very well in the central region (\y{\ < 1.1), which is the region used in the W 
mass determination 0. The difference between the two calculations is mainly at \y{\ > 1.7, 
as illustrated in Figure [p. 

In summary, the asymmetry data demonstrate the value of collider data in the measure- 
ment of parton distribution functions and place the strongest constraint on the - ratio of 
quark momentum distributions in the proton over the range x of 0.006 to 0.34 at Q 2 « Mjy. 
The data indicate the need for modifications in the - ratio in the PDFs constrained by the 
previous CDF asymmetry data. 
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TABLES 



TABLE I. The mean x values of the d and u parton distributions probed by the CDF asymmetry 
data for the W + boson as a function of the lepton rapidity. The | ratios from the CTEQ-3M 
parameterization are shown for these x values. The mean x values for the W~ at positive y\ are 
the same as those for the W + at the same negative y\. 



ton) 


{x u ) 




iM 




-2.20 


0.009 


0.219 


0.911 


0.412 


-1.54 


0.016 


0.129 


0.864 


0.542 


-1.10 


0.021 


0.088 


0.837 


0.625 


-0.49 


0.038 


0.051 


0.765 


0.722 


n nn 
u.uu 


U.U44 


n fi/i/i 

U.U44 


0.744 


0.744 


U.4y 


n non 
u.uyy 


n mi 
U.Uzl 


0.601 


0.837 


1 in 
1.1U 


n ice 


U.Ulz 


0.483 


0.889 


1 KA 
1.04 




n nnn 

u.uuy 


0.405 


0.911 


i on 


U.OOO 


U.UUO 


0.306 


0.936 


TABLE II. 


Backgrounds (%) 


in the W —> eu and W - 


-> /ii/ charge asymmetry event samples. 


Sample 




QCD dijets 


Z 


e + e or fi + fi 


Central e 




0.7 ±0.2 




< 0.2 


Central [i 




0.6 ±0.2 




4.7 ±0.7 


Plug-SVX e 




1.6 ±0.3 




< 0.2 


Plug-CTC e 




2.4 ±0.6 




< 0.2 


Forward \i 




4± 2 




6.5 ±0.6 
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TABLE III. The charge asymmetries (after all corrections) and statistical and systematic un- 
certainties in the combined e and fx channels. 



(\yi\) 


A(yi) 


°"sys 


^stat+sys 


0.11 


0.017 


±0.000 


±0.008 


0.30 


0.050 


±0.000 


±0.007 


0.50 


0.097 


±0.000 


±0.007 


U. 1 u 


U. lOZ 


_i_n nm 


_i_n nns 

ItU.UUo 


0.89 


0.138 


±0.001 


±0.009 


1.05 


0.163 


±0.002 


±0.016 


1.31 


0.122 


±0.002 


±0.013 


1.54 


0.121 


±0.003 


±0.012 


1.79 


0.075 


±0.005 


±0.023 


1.95 


0.025 


±0.010 


±0.029 


2.20 


-0.095 


±0.015 


±0.023 
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FIGURES 
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FIG. 1. The fully corrected charge asymmetry. Data from all the detectors for positive and 
negative yi are combined. The statistical and systematic errors are added in quadrature. 
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